Phase Transition Properties of 3D Potts Models 
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Using multicanonical Metropolis simulations we estimate phase transition prop- 
erties of 3D Potts models for q = 4 to 10: The transition temperatures, latent 
heats, entropy gaps, normalized entropies at the disordered and ordered endpoints, 
interfacial tensions, and spinodal endpoints. 

I. INTRODUCTION 
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Potts models were introduced as a footnote in the thesis by Potts [1| on clock models, 
whose investigation had been proposed to him by his adviser Domb. We define their energy 
■ function by 

Si £ (fe) = 2£(i-%f\gf)) (1) 

; (ij) \ q > 

^ ! where 5(qi,qj) is the Kronecker delta function, the sum (ij) is over the nearest neighbor 
lattice sites, the index (k) refers to the configuration, and gj is the Potts spin at site i. 
For the g-state model Potts spins take the values 0, 1, . . . , q — 1. Following the notation of 
Ref. [2j the factor of two in front of the sum is introduced to match for q = 2 on Ising model 
conventions. Expectation values of the Gibbs canonical ensemble are calculated with the 
Boltzmann (or Gibbs) factor exp(— /3E), (3 = l/(kT). 

Potts, and not so much clock models, received considerable attention up to the day. 
Developments till 1982 were reviewed by Wu {^J. In 2D the models have second order phase 
transitions for q = 2 to 4, and first order transitions for q > 5. In a work by Baxter 0], 
critical temperatures, latent heats and entropies were analytically derived. Equations for 
interface tensions followed later [5[ and played a role in the verification of the multicanonical 
approach \^ to Markov Chain Monte Carlo (MCMC) simulations. In 3D the q = 2 Ising 
model phase transition is second and the q > 3 transitions are first order. The strength 
of the first order transitions, measured by their latent heat, increases with q. The q = 5 
transition in 2D and the q = 3 transition in 3D are weakly first order. 

Potts models proved amazingly versatile to grasp the essence of physically interesting 
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situations, many of them are described in Ref. [3j. When they are generalized by introducing 
quenched random variables as exchange constants, the Ising case of the Edwards-Anderson 



spin glass [7j and quadrupolar Potts glasses [8] are obtained. Particular choices give anti- 
ferromagnets and the fully frustrated Ising model 0]. 

The 3D q = 3 Potts model shares the center symmetry of SU(3) gauge theory [10'] and 
plays a role in our understanding of the QCD deconfining phase transition. This and other 
motivations led to a number of numerical investigations of the 3D 3-state Potts model 
[ill . [12I I3I Q, Ifil [l7, [l8]. However, when we recently looked out for a strong first 
order transition in 3D, we found only a few papers 19|, |2Qi, |2l| dealing with 3D Potts for 
q > 4. The purpose of this paper is to fill this gap in the literature for 3D Potts models 
up to q = 10. Thereby, we will closely follow the outline of a previous investigation of 
the equilibrium statistical properties of the 3D 3-state Potts model by two of the present 
authors (l8| . 

Next, we briefly summarize our simulation method and give an overview of our assembled 
statistics. In section [TTT1 we calculate and analyze transition temperatures, specific and latent 
heats. Section IIVI deals with entropy and energy across the phase transition, section |V] with 
interface tensions and spinodal endpoints. A brief summary and conclusions are given in 
the final section I VI I 



II. SIMULATION METHOD AND OVERVIEW OF DATA 

We want to calculate expectation values O in the Gibbs canonical ensemble. They are 
given by the ensemble average 

K 

6 = d(f3) = (0)f 3 = Z- l J20 {k) e- /3E(k) (2) 

k=l 

where 

Z = Z{f3) = jZe-? E{k) =Y.n{E)e-? E (3) 

k=l E 

is the partition function. The index (k), k = 1, . . . , K labels the configurations (microstates) 
of the system and E^ k ' is the (internal) energy of configuration {k). In the last equation 
n(E) is the number of configurations with energy E. 

We consider Potts models on cubic lattices of size L 3 with periodic boundary conditions. 
There are N = L 3 Potts spins. Each microstate k defines a unique arrangement of Potts 
spins and vice versa: 

k = {q{ k \...,qP}. (4) 
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As each Potts spin can take on q values, there are 

K L = Z L (0)=Y,n L (E) = q N , (5) 

E 

microstates. Even for rather small numbers of L, K is a very large, so that one will not be 
able to sum the partition function explicitly. Instead, one can use statistical methods. 

MCMC simulations are a suitable approach to evaluate equilibrium properties of the 
canonical ensemble. Off the phase transition temperatures canonical Metropolis or heatbath 
simulations with weight factor exp(— f3E) work sufficiently well, provided a ordered start is 
used for simulations in the ordered phase, and (somewhat less important) a disordered start 
for simulations in the disordered phase. However, canonical simulations deteriorate quickly 
when it comes to the (most interesting) investigation of phase transition properties. For 
first order transitions the reason is that the relevant transition states are in the canonical 
ensemble suppressed ~ exp(— 2ondL D ~ 1 ), where a d is the order-disorder interface tension. 

Multicanonical simulations [2, (J are a remedy for this supercritical slowing down. One 
performs for E m \ n < E < E max MCMC simulations with a working approximation of the 
weight factors 

w m UE)=e- b ^ E+a ^ = ^-, (6) 

n(E) 

supplemented by the canonical weights exp(— P ma _ x E) for E < E min and exp(— /3 min £') for 
E > E max . Here n(E) is the number of states with energy E as introduced by Eq. Q. 
With the weights ([6]) all energies in the range E m i n < E < E max are sampled with the same 
probability, so that the Markov process will perform some kind of random walk in the range 
-E'min < E < £? mH . The microcanonical inverse temperature b(E) and the dimensionless free 
energy a(E) follow uniquely from n(E) due to the relation 

a(E -e) = a(E) + [b(E - e) - b(E)] E , a(£ max ) = , (7) 

where e is the step to the next energy. However, n(E) is a-priori unknown, so that some 
iterative procedure needs to be used. In this context working approximation of the weights 
means that any weights that enable cycling (also called tunneling) 

£ max -> ^min and back (8) 

are considered to be acceptable. Actually it is known that not the weights but some 
modifications of them, can be optimal for that purpose 
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there is a residual exponential slowing down 241 ] . 



23| and it should be noted that 



We find working estimates of the multicanonical weights by finite size (FS) extrapolations 
from a smaller lattice to a next larger lattice, a method which was already used in 0]. 
Obviously, this requires that the FS behavior of the system under consideration is well 
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defined. In particular, for complex systems like spin glasses orproteins this is not the case 
and more sophisticated recursion approaches need to be used [2J , most noted is presumably 



the one by Wang and Landau [25]. 



Let e = E/N = E/L 3 be the energy density. For our Potts models we simply convert the 
microcanonical inverse temperature b^E) for a given lattice size L\ to the b^{E) of our 
next larger lattice size L2 > L\ via the interpolation 

b L2 (E) = w + b Ll (E+) + u>_ b Lx (£_) with w ± = \ 6 ~ 6 ^ , (9) 

6_|_ — 6_ 

where e = _E/(L 2 ) 3 and e_ < e and e+ > e are the values closest to e, so that we have for the 
corresponding E + and E_ values entries in b^E). For our purposes this simple procedure 
turned out to be sufficiently accurate. Better extrapolations can be expected by taking 
into account details of the shape of til{E), as for a magnetic field driven phase transition 



discussed in Ref. [26] 



Suitably, the E max , E min , /3 min and /3 max parameters, which accompany the weights (jUJ) 
are chosen so that they embrace the phase transitions and that 

-Emax = (E)(3 mia , E min = (-E)/3 max (10) 

holds. From a simulation with these weights, canonical expectation values are obtained by 
reweighting for the temperature range /3 min < (3 < /5 max . This property has coined the name 
multicanonical. To calculate the partition function ([3]), from which the normalized entropy 
and free energy follow, one has to include (3 — in this temperature range. So we choose 
Anin = 0, for which our normalization of the energy (CQ) implies E mSuX = 0. For j3 max we chose 
the values given in table [U each of them well above the transition value (3 C . Our simulations 
of multicanonical ensembles defined by the weights (EI) rely on the Metropolis algorithm. 
We update sequentially with one Metropolis update per spin during one sweep through the 
lattice. This is more efficient 0] than picking spins at random for the updates. 

In table fl] we give an overview of the statistics per run. We followed the outline of Potts 
model MCMC simulations in 0]. First, we performed the number of sweeps listed in table [I] 
for reaching equilibrium. Data from these sweeps are excluded from the statistics for which 
measurements were performed. Subsequently, we collected for each run 32 histograms, each 
relying on the number of sweeps listed in the table. All error bars are then calculated 
with respect to these 32 bins (32 jackknife bins when nonlinear operations on the data are 
involved). From the student distribution it is known that error bars from 32 independent 
Gaussian data give almost Gaussian confidence probabilities at the level of two standard 
deviations. 

For most data points we performed two runs with the statistics of table [B The first runs 
are based on the weights iterated from the closest smaller lattice. These data are taken 
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TABLE I: Analyzed statistics per production run in sweeps: 32 x the number given. 



L 


q = A 


q = 5 


q = 6 


q = 7 


q = 8 


q = 9 


q = 10 


/^max — 


0.45 


0.41 


0.43 


0.45 


0.47 


0.49 


0.50 


2 


256 


1024 


6xl0 3 


10 x10 s 


15xl0 3 


15xl0 3 


16384 


3 


1024 


2048 


9x 10 3 


10 x10 s 


15x 10 3 


15xl0 3 


4000 


4 


1024 


2048 


9xl0 3 


10x10 s 


20xl0 3 


20xl0 3 


4000 


6 


4096 


8192 


lxlO 4 


15xl0 3 


20xl0 3 


30xl0 3 


65536 


8 


4096 


32768 


4xl0 4 


45xl0 3 


5xl0 4 


75xl0 3 


262144 


10 


16384 


lOxlO 4 


16xl0 4 


25xl0 4 


3xl0 5 


4xl0 5 


1048576 


12 


65536 


25xl0 4 


3xl0 5 


45xl0 4 


55xl0 4 


7xl0 5 


4194304 


14 


262144 


45xl0 4 


5xl0 5 


14xl0 5 


17xl0 5 


22xl0 5 


16777216 


16 


262144 


85xl0 4 


15xl0 5 


20xl0 5 


26xl0 5 


32xl0 5 


2xl0 7 


18 


1048576 


13xl0 5 


21xl0 5 


29xl0 5 


70xl0 5 


2xl0 7 


4x10' 


20 


4194304 


19xl0 5 


30xl0 5 


50xl0 5 


14xl0 6 


8xl0 7 


8xl0 7 


99 
zz 








1 Y in7 

1 A 1U 


O A 1U 






24 


45xl0 5 


39xl0 5 


lxlO 7 


3xl0 7 








26 


9xl0 6 




2xl0 7 










28 


14xl0 6 


lxlO 7 












30 


15xl0 6 


3xl0 7 













to refine the weights for the lattices at hand. The refined weights are used for the second 
production runs on these lattices. Exceptions from this procedure are iterations from a 
smaller to a larger lattice immediately after the first run. This speeds up the process of 
getting to larger lattices and has often been done when the cycling frequency of the first run 
was already satisfactory. 

Table [IT] collects the number of cycling (jSJ) events obtained in the first production runs. 
For L = 2, some cycling is already achieved by a canonical simulation at (3 = 0. This allows 
one to determine multicanonical weights for the second run on 2 3 lattices, which have large 
cycling rates as shown in table IHIl and to start off FS iterations of the weights. From the 
2 3 lattices we extrapolate weights for the first runs on 3 3 lattices, refine them for the second 
runs on 3 3 lattices, iterate to the next larger lattice, and so on (up to the before mentioned 
exceptions). The calculations were carried out on PC clusters at FSU. Our present lattice 
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TABLE II: Number of cycling events for the first production runs. 





g = 4 


q = 5 


q = 6 


q = 7 


q = 8 


q = 9 


q = 10 


L = 2 : 


3 


6 


4 


2 


1 


1 


1 


L = 3 : 


32 


135 


350 


246 


239 


104 


15 


L = 4 : 


109 


139 


377 


358 


550 


531 


88 


L = 6 : 


68 


117 


80 


98 


89 


87 


164 


L = 8 : 


21 


125 


106 


76 


71 


92 


113 


L = 10 : 


20 


186 


190 


189 


122 


109 


135 


L = 12 : 


77 


249 


195 


169 


99 


52 


107 


L = 14 : 


173 


279 


178 


242 


102 


42 


86 


L = 16 : 


125 


313 


301 


139 


47 


12 


11 


L = 18 : 


353 


322 


228 


69 


38 


10 


3 


L = 20 : 


711 


327 


172 


48 


11 


5 


2 


r 99 . 

X/ — ZZ . 








97 


D 






L = 24 : 


441 


84 


1 


27 








L = 26 : 


903 




116 










L = 28 : 


867 


43 












L = 30 : 


356 


639 













sizes are limited by the deterioration of cycling with increasing L, the computational power 
of a single PC (here 2-3 GHz per PC), and the limitation of the total length of one run to 
a few months. 

During the simulations we collect histograms of the energy in the multicanonical ensem- 
ble, h mu (E), and calculate functions of the energy f(E) from them by reweighting to the 
canonical ensemble: 

= Ee /Cg) Ug e*p [-13 E + b(E) E - a(E)] 
1 J2 E hmu(E) ex V [-PE + b(E)E-a(E)} ' 1 ' 

where the sums are over all energy values for which h mu (E) has entries. The computer 
implementation of this equation requires care, because the differences between the largest 
and the smallest numbers encountered in the exponents can be large. We rely here on the 
logarithmic coding of Ref. [2]. Whenever the function f(E) is non-linear jackknife binning 
is employed. 
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TABLE III: Number of cycling events for the second production runs. 





g = 4 


q = 5 


q = 6 


q = 7 


q = 8 


q = 9 


q = 10 


L = 2 : 


139 


1451 


5918 


7091 


2614 


7985 


6339 


L = 3 : 


288 


378 


1275 


1167 


1436 


1209 


1098 


L = 4 : 


72 


126 


410 


331 


578 


497 


324 


L = 6 : 


69 


100 


109 


127 


117 


165 


311 


L = 8 : 


22 


154 


125 


117 


105 


110 


349 


L = 10 : 


42 


180 


223 


295 


243 


273 


598 


L = 12 : 


93 


256 


221 


228 


218 


246 


1142 


L = 14 : 


209 


251 


219 


427 


384 


394 


2423 


L = 16 : 


130 


308 


366 


396 


352 


359 


1542 


L = 18 : 


323 


301 


347 


342 


585 


1056 


1707 


L = 20 : 


921 


307 


329 


358 


590 


2168 


1591 


r 99 . 

X/ — ZZ . 










770 
(It) 






L = 24 : 


521 


321 


487 


910 








L = 26 : 






708 










L = 28 : 


986 


460 












L = 30 : 


869 















III. ENERGIES, TRANSITION TEMPERATURES AND LATENT HEATS 

In Fig. [1] we show internal energies as functions of j3, for each q from the largest available 
lattice (error bars are not resolved on the scale of this figure) . This gives a rough estimate of 
the latent heats and the inverse transition temperatures (3 t . Accurate results follow from FS 
extrapolations of indicators, which are defined on finite lattices, so that they converge in the 
limit L — > oo (quickly) towards the infinite volume value of the desired physical quantity. 

We calculate specific heats via the fluctuation-dissipation theorem 

C = {(E 2 ) - (Ef) . (12) 

For first order phase transitions the finite volume specific heats are regularization of Dirac 
delta functions, which are the infinite volume extrapolations. The multicanonical approach 
allows to calculate specific heat values for a continuous range of (3 values. Consequently, the 
locations of the maxima can be accurately determined. 



8 



LU 







-1 



-5 




2 


- q= 


=10, L= 


=20 




q= 


=09, L= 


=20 




q= 


=08, L= 


=22 


3 


q= 


=07, L= 


=24 




q= 


=06, L= 


=26 




q= 


=05, L= 


=30 


4 


" q= 


=04, L= 


=30 




q= 


=03, L= 


=50 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
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FIG. 1: Internal energies. 



Finite lattice indicators for the transition temperatures are called pseudo transition tem- 
peratures, and there are various options to define them. On finite volumes their values differ, 
while they all converge to the same L — > oo limit. We use here three definitions of pseudo 
transition temperatures: /5£ t , the j3 value at which equal heights are achieved in the double 
peaked energy histogram, (3p t) the position of the central energy of the latent heat, and (3p t , 
the location of the maximum of the specific heat. The first two definitions are explained 
in more detail later, (3^ t plays a role in determining the entropy gaps in section [TV] and /3* t 
for the interface tensions in section |V] (they are labeled in this order for consistency with 
Ref. Q). 

For q = 4 Fig. [2] shows our linear fits 



fU L ) = % + 



L 3 



(13) 



which determine the infinite volume transition temperatures. The vertical order of the 
fits agrees with that in the legend. While the finite volume estimators (3 l pt {L) differ, the 
infinite volume extrapolations are consistent with one another. The smallest lattices have 
been omitted from the fits to ensure an acceptable goodness of fit Q in each case. This 
q = 4 pattern repeats for all q: The quality of the fits from our three definitions of pseudo 
transitions temperatures are similar, and the final estimates as well as their error bars are 
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CD. 



0.3144 




0.3142 




0.314 




0.3138 




0.3136 


- P 


0.3134 


- P 


0.3132 


- P 


0.313 






Ik 



0.0002 0.0004 0.0006 0.0008 0.001 

L- 3 

FIG. 2: Transition temperature fits for q = 4. 



consistent with one another. To give one best number for each q, we simply average over the 
three estimates. We average also their error bars, because all three estimators rely on the 
same simulation, so that one does not expect error bar reduction when averaging over them. 



e IIVI For the convenience of 
li and 



27| . respectively So 



The thus obtained transition values (3 t (q) are collected in tab 
the reader we have included q = 3 and 2 estimates from Ref. 
far we have not found a simple formula for the q dependence like fit = ln(l + v /g)/2, which 
holds in 2D |4J (given here in our convention). 

Following Ref. 28(], we extract the latent heat Ae = AE/N by fitting the maxima of the 
specific heat to the form 

C max (L) = ai + a 2 L 3 (14) 

and using the relation 

a 2 = i(A) 2 (Ae) 2 . (15) 



The results are also included in table IIVI 



IV. ENTROPY AND ENERGY ACROSS THE PHASE TRANSITIONS 



The entropy density is 



(16) 
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TABLE IV: Estimates of observables (q = 2 from Ref. [27J] and q = 3 from Ref. 181 ]) 



Q 


A 


Ae 


e(A) 


e+ 


e 


2 


0.2216544 (06) 





-0.9957 (14) 


e(A) 


e(Pt) 


3 


0.2752827 (29) 


0.3286 (15) 


-1.3470 (74) 


-1.1826 (73) 


-1.5112 (79) 


4 


0.3143041 (17) 


1.16294 (61) 


-1.719 (36) 


-1.1367 (64) 


-2.3019 (64) 


5 


0.3447205 (12) 


1.84619 (20) 


-1.987 (20) 


-1.063 (20) 


-2.910 (20) 





U.369lWU (15J 


z.3644z (1/J 


-2.1 ( / (26 J 


n one; / 
-U.y95 (26J 


-3.35y (26J 


7 


0.3909657 (17) 


2.76430 (12) 


-2.316 (22) 


-0.934 (22) 


-3.698 (22) 


8 


0.4094959 (23) 


3.08039 (15) 


-2.421 (32) 


-0.881 (32) 


-3.961 (32) 


9 


0.4259432 (23) 


3.33628 (12) 


-2.503 (32) 


-0.835 (32) 


-4.171(32) 


10 


0.4407371 (18) 


3.547570 (87) 


-2.567 (25) 


-0.794 (25) 


-4.341 (25) 



where / is the free energy density, which is continuous at the phase transition. So, the 
entropy gaps across the phase transitions are 



As = fa Ae. 



(17) 



or As = 2y / 7i^ with a,2 from the fits ( tHj) . The entropy and energy density endpoints in the 
disordered (+) and ordered (— ) phases are given by 



(3t(e + -f([3 t )) 
A(e--/(ft)) 



(18) 
(19) 



and are more difficult to compute than the gaps, because the additive normalization con- 
stants no longer drop out. We follow the method of Ref. [18|, which relies on the definition 
of Pp t as given below. 

In multicanonical simulations, the normalization constant for the entropy is determined 
by the known value at /3 = 0: 



On 



S = In (q N ^J and s = = Inq . 



(20) 



Fig. [3] shows normalized entropy densities of our models for our largest lattices. 

To calculate the endpoints of the entropy and energy on the ordered and disordered sides 
of the transitions, we define fip t (L) by the relation 

1 

2 



e L (f3 2 pt ) 



eiifc) + e£(J3. 



(21) 



CO 



2.5 



1.5 



0.5 







ll 



q=10, L=20 
q=09, L=20 
"q=08, L=22 
q=07, L=24 
q=06, L=26 
"q=05, L=30 
q=04, L=30 
q=03, L=50 







0.1 



0.2 



0.3 



0.4 



0.5 



P 



FIG. 3: Normalized entropy densities. 



LU 
Q_ 




-0.5 



FIG. 4: Energy density histograms at (3p t {L) for q = 4. 
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TABLE V: Estimates of observables (q = 3 from Ref. [181] ) . 



q 




f(fit) 




s 


w + 


w 


2 


0.55715 (31) 


-3.50956 (23) 


s(A) 


a(A) 


1.75 


1.75 


3 


0.8491 (21) 


-4.431364 (50) 


0.8943 (21) 


0.8038 (22) 


2.45 


2.24 


4 


0.983 (11) 


-4.846358 (41) 


1.166 (11) 


0.800 (11) 


3.21 


2.23 


5 


1.0680 (70) 


-5.084679 (34) 


1.3862 (70) 


0.7498 (70) 


4.00 


2.12 





1.1324 (93 J 


-5.239933 (50J 


1.0D95 (93 J 


n fin A ( no\ 
0.6954 (93) 


A on 
4.80 


2.00 


7 


1.1858 (85) 


-5.349346 (60) 


1.7262 (85) 


0.6455 (85) 


5.62 


1.91 


8 


1.233 (13) 


-5.430753 (66) 


1.863 (13) 


0.602 (13) 


6.44 


1.83 


9 


1.274 (14) 


-5.493609 (42) 


1.984 (14) 


0.563 (14) 


7.27 


1.76 


10 


1.312 (11) 


-5.543856 (50) 


2.094 (11) 


0.530 (11) 


8.11 


1.70 



where e^((3p t ) are the locations of the maxima of the double peak histogram at f3^ t {L). For 
q = 4 these histograms are shown in Fig. [4] (we excluded the L = 30 lattice to keep a 
reasonable scale in the figure). This construction ensures that the energy endpoints ef are 
positioned symmetrically about the central energy density ei(/3p t ): 



4 = e L (fa) ± \/\e L (22) 



and one finds that 



holds as well. As in 



4 = s L ((3 2 pt ) ± l -/\s L (23) 

isl we use jackknife estimators and arrive at the values for e(/5j), e + and 
compiled in table IIVI and those for s(fi t ), fiftt), s + and s~ compiled in table [V] Although 
there are simple relations between these values, we have to list them separately, because 
these relations do not determine error bars. The estimates are correlated and the jackknife 
procedure takes care of correct error bars. Also it should be noted that (3 t in the arguments 
of e, s and / is /? t 2 , defined as the extrapolation of the pseudo transition temperatures as 
defined by fj21jl . These values are consistent with the (3t values listed in table HVt but not 
identical, as fit of table IIVI is the average of the extrapolations from our three definitions of 
pseudo transition temperatures. 

In the last two columns of table |V] we give the first few digits of the effective number of 
states on the disordered and ordered sides of the phase transition, 

w ± = exp(s ± ) . (24) 
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FIG. 5: Equal heights energy density histograms q = 10. 



Amazingly, the effective number of states per spin at the the ordered endpoint goes down 
by increasing q. 



V. INTERFACE TENSIONS AND SPINODAL ENDPOINTS 



For L — > oo the interface tension between ordered and disordered phases is 29] 

2a od (L) = ^ ^ f^H^l (25) 



where P max (L) represents the value of the maxima when the energy histogram is reweighted 
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to equal heights and P m j n (L) the minimum in between the peaks, 
equal heights histograms in Fig. [51 Including capillary waves 
3-parameter fits to the form (compare Eq. (16) of [33]) 

2 a od (L) + -^j- = 2a od + — + — 



7 or q = 10 we show our 
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321 ] , we perform 2- and 



(26) 



In case of the 2-parameter fits we set c 3 = 0. While the 3-parameter fits are somewhat 
unstable, consistent 2-parameter fits are limited to the largest three lattices. The differences 
between these fits exhibit systematic errors, which show that larger lattices are needed for 
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TABLE VI: Estimates of interface tensions (yes/no refers to capillary waves, q = 3 from Ref. [18I]). 



Q 


2 o- orf - yes 


2 a a d - no 


Q 


2 cr od - yes 


2 a d - no 


3 


0.001806 (35) 


0.001602 (35) 


7 


0.1484 (15) 


0.1478 (17) 


4 


0.0224 (11) 


0.0221 (13) 


8 


0.1897 (56) 


0.1891 (56) 


5 


0.0632 (21) 


0.0628 (22) 


9 


0.2308 (87) 


0.2302 (84) 


6 


0.1054 (50) 


0.1050 (45) 


10 


0.2628 (40) 


0.2688 (47) 



0.48 
0.47 
0.46 
0.45 
cQ- 0.44 
0.43 
0.42 
0.41 
0.4 

FIG. 6: Microcanonical temperature b(E) for q = 10. 




high precision results. The results of the 2-parameter fits are compiled in table IVIl where 
the differences to the less stable 3-parameter fits are used to estimate systematic errors, 
which are, in these cases, substantially larger than statistical errors of the fits. 

We also include in table IVIl results from fits without the capillary wave contribution 
ln(L)/(2L 2 ). With the exception of the q = 3 case from Ref. [18] , the difference between the 
two fits is always smaller than the expected error from other sources. 

The double peak histo gram s at first order phase transitions are intimately related to a 



Maxwell construction [34J, [35] for the inverse microcanonical temperature b(E) defined by 
Eq. (jSj). For q — 10 and two lattice sizes this is shown in Fig. [61 The areas above and below 
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the infinite volume line are identical and for L = 20 one sees that a small fraction of the 
curve joins this line. 

The minimum (3^ m and the maximum (3^ &x of the (3(E) curve are the inverse spinodal 
temperatures. Equilibration at (3 with dissipative model A (Glauber) dynamics [36| encoun- 
ters metastability in the range (3 t < (3 < (3^ x after a disordered start, whereas after an 
ordered start it encounters metastability in the range (3 t > (3 > P^m- F° r the (3(H) of a 
magnetic field driven phase transition 26J this would already be the entire metastability 
picture. In case of the temperature driven phase transitions of Potts model it is more com- 
plicated, because metastability after a disordered start persists for equilibration at f3 > P^ ax 
due to order-order domain walls, which are for q = 3 investigated in Ref. 37]. 

In the past there may have been some hesitation in identifying (3^ in and (3^ ax , as defined 
here, with the spinodal endpoints. The reason is that their values agree in the infinite volume 



limit with (3 t |34|, |35|] as is illustrated by the dotted line in Fig. El So the metastability 



disappears in the infinite volume limit, whereas the opposite is the case for the mean field 
spinodal, which is introduced in many textbooks 38[]. However, the recent finite volume 
analysis [39] of Kolmogorov-Johnson-Mehl-Avrami (KJMA) theory demonstrates that the 
mean field approach is a conceptually wrong starting point for describing the infinite volume 
limit of phase conversion. Within the KJMA framework one gets for V — » oo always spinodal 
decomposition [40J and never metastability. Our definitions of (3^ iri and (3^ &x are consistent 
with this picture as well as with studies of magnetic field driven phase transitions by Rikvold 
et al. [4l|. 

For large L the areas in the Maxwell construction are known 35| to shrink ~ 1/L. 
Therefore, the leading order 2-parameter fit for (3 sp (L) is 



(3 sp (L)-(3 t 



(27) 



Using (3 t from table HVl we show in Fig. [7] the fits to this form. Together with their goodness 
of fit Q 2] the fit parameters are collected in table IVHI The Q values are a bit on the 
high side, as a relatively flat (3(E) curve tends to give rather large statistical errors for the 
spinodal estimates (3 sp (L). In table IVHl this is reflected by 02 parameters, which are mainly 



statistical noise about zero. For the q = 3 data [18j, (3(E) is altogether too flat to allow for 



reasonably accurate (3 sp (L) estimates (larger lattices would be needed). 



VI. SUMMARY AND CONCLUSIONS 



For 3D, g-state Potts models in the range q = 4, . . . , 10 we have estimated a number 
of observables by multicanonical MCMC calculations and supplemented them with q = 2 
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FIG. 7: FS behavior of the spinodals. 



TABLE VII: Estimates of the fit parameters of Eq. (|27|h 



q 




Q>2,max 


Q 




^2,min 


Q 


4 


0.0193 (21) 


1.6 (3.0) 


0.57 


-0.0053 (24) 


21 (21) 


0.83 


5 


0.0357 (20) 


3.0 (1.5) 


0.19 


-0.0154 (19) 


5.5 (3.7) 


0.84 


6 


0.0543 (20) 


2.07 (89) 


0.84 


-0.0200 (29) 


7.5 (4.3) 


0.81 


7 


0.0764 (22) 


0.53 (57) 


0.67 


-0.0315 (27) 


2.3 (2.0) 


0.95 


8 


0.0969 (10) 


-0.01 (19) 


0.85 


-0.0412 (27) 


0.8 (1.4) 


0.69 


9 


0.1154 (19) 


-0.14 (30) 


0.58 


-0.0475 (25) 


0.6 (1.1) 


0.82 


10 


0.13631 (93) 


-0.58 (12) 


0.05 


-0.0639 (19) 


-2.00 (41) 


0.21 



and q = 3 results from the literature. Transition temperatures, latent heats and energy 
endpoints of the phases are given in table IIVt entropy and free energy values in table |Vj 

Less accurate are our interface tension estimates of table |VT1 They could possibly be im- 
proved by using simulation techniques similar to those, which led to high-precision estimates 
of the order-order interface tension in the 3D Ising model 42]. 

Minima and maxima of the microcanonical inverse temperature curve b(E) are identified 
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as adequate definition of spinodal endpoints. As expected [39j, |4l| the thus defined regions 
of metastability disappear in the infinite volume limit. 

We hope that future investigations of 3D first order phase transitions will benefit from 
the results collected in this paper. 
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